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ABSTRACT 

Massive young clusters (YCs) are expected to host intermediate-mass black 
holes (IMBHs) born via runaway collapse. These IMBHs are likely in binaries and 
can undergo mergers with other compact objects, such as stellar mass black holes 
(BHs) and neutron stars (NSs). We derive the frequency of such mergers starting 
from information available in the Local Universe; in particular, we assume that 
a fraction ~ 0.75 of all the YCs more massive than 10 5 M might host one 
IMBH, as suggested by a statistical analysis of the properties of YCs in the 
Milky Way and in the Antennae. Mergers of IMBH— NS and IMBH— BH binaries 
are sources of gravitational waves (GWs), which might allow us to reveal the 
presence of IMBHs. We thus examine their detectability by current and future 
GW observatories, both ground- and space-based. In particular, as representative 
of different classes of instruments we consider Initial and Advanced LIGO, the 
Einstein gravitational-wave Telescope (ET) and the Laser Interferometer Space 
Antenna (LISA). We find that IMBH mergers are unlikely to be detected with 
instruments operating at the current sensitivity (Initial LIGO). LISA detections 
are disfavored by the mass range of IMBH— NS and IMBH— BH binaries: less than 
one event per year is expected to be observed by such instrument. Advanced 
LIGO is expected to observe a few merger events involving IMBH binaries in 
a 1-year long observation. Advanced LIGO is particularly suited for mergers of 
relatively light IMBHs (~ 10 2 M ) with stellar mass BHs. The number of mergers 
detectable with ET is much larger: tens (hundreds) of IMBH-NS (IMBH-BH) 
mergers might be observed per year, according to the runaway collapse scenario 
for the formation of IMBHs. We note that our results are affected by large 
uncertainties, produced by poor observational constraints on many of the physical 
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processes involved in this study, such as the evolution of the YC density with 
redshift. 

Subject headings: gravitational waves - - black hole physics - - galaxies: star 
clusters: general — (stars:) binaries: general — stars: kinematics and dynamics 

1. Introduction 

So far, there are no definitive observational proofs of the existence of intermediate-mass 
black holes (IMBHs), i.e. of black holes (BHs) with a mass ranging from ~ 10 2 to ~ 10 5 M 
(see van der Marel 2004 for a review). However, the formation of IMBHs is predicted by 
various theoretical scenarios. In the early Universe, IMBHs might form via the direct collapse 
of metal-free stars (Heger et al. 2003) or via the collapse of gaseous disks in the center of pre- 
galactic halos (Begelman, Volonteri & Rees 2006). IMBHs might form even in more recent 
epochs, via the direct collapse of metal-poor stars (Mapelli, Colpi & Zampieri 2009; Zampieri 
&: Roberts 2009; Mapelli et al. 2010), the repeated mergers of stellar mass BHs in globular 
clusters (GCs, Miller & Hamilton 2002) and the runaway growth of IMBHs (Portegies Zwart 
& McMillan 2002, hereafter PZM02) in massive (> 1O 4 M ) young (< 3 x 10 7 yr) clusters 
(YCs). In particular, the runaway growth scenario predicts that star clusters with initial half- 
mass relaxation times < 25 Myr are dominated, during the core collapse, by star collisions, 
which build up one or even two (Giirkan, Fregeau & Rasio 2006) very massive objects, which 
likely evolve into IMBHs (PZM02). Thus, YCs are among the best environments where we 
can search for IMBHs. 

However, the observational features of IMBHs are difficult to pinpoint. There are hints 
(Strohmayer & Mushotzky 2003; Kaaret, Ward & Zezas 2004; Strohmayer et al. 2007) 
that IMBHs might power some of the brightest ultra-luminous X-ray sources (ULXs, i.e. 
those sources with a X-ray luminosity higher than the one predicted for a ~ 10 M Q BH). 
Other studies point out that the presence of IMBHs may be inferred from anomalies in the 
acceleration of millisecond pulsars (Colpi, Possenti & Gualandris 2002; Colpi, Mapelli & 
Possenti 2003) or from the observation of high- velocity stars in GCs (Mapelli et al. 2005). 
Furthermore, the combination of kinematic and photometric data with dynamical models 
suggests the presence of IMBHs in the GC Gl (Gebhardt, Rich & Ho 2002, 2005) in the 
Andromeda galaxy and in the GC M15 (Gerssen et al. 2002; van den Bosch et al. 2006) in 
the Milky Way (MW). 

Finally, gravitational waves (GWs) will provide a powerful tool to detect IMBHs. In 
fact, IMBHs lying in the center of dense clusters are expected to form binaries with stars or 



-3 - 



stellar mass BHs (e.g. Sigurdsson & Hernquist 1993). Such binaries harden (i.e. progressively 
reduce their orbital separation) due to gravitational encounters (Colpi et al. 2003 and 
references therein). When the orbital separation is sufficiently small and if the companion 
of the IMBH is a compact object, such as a stellar mass BH or a neutron star (NS), the 
binary enters a regime in which the orbital evolution is driven by radiation reaction induced 
by GW emission. This further reduces the orbital separation, until the system merges. GWs 
emitted by binaries hosting IMBHs are mostly in the frequency range accessible to current 
or forthcoming detectors, such as the Laser Interferometer Gravitational- Wave Observatory 
(LIGO; Barish & Weiss 1999; Abbott et al. 2009) and Virgo (Acernese et al. 2004) in initial 
and advanced configuration (also known as first and second generation laser interferometers). 
GWs from IMBH mergers also fall in the frequency range of future GW observatories, such 
as the Laser Interferometer Space Antenna (LISA; Bender et al. 1998) and third generation 
ground-based instruments, for which we adopt, as an example, the concept of the Einstein 
gravitational- wave Telescope (ET; Freise et al 2008, Hild et al. 2008). In particular, Miller 
(2002), considering IMBHs in GCs, predicts that tens of sources connected with Galactic 
IMBHs will be detectable with LISA in a 5 yr integration. More recent studies investigate 
the emission of GWs from IMBH— IMBH binaries (Fregeau et al. 2006), whose existence is 
predicted by simulations (Giirkan et al. 2006). Massive stellar mass BHs (15 — 20 M ) are 
also found to be important sources of GWs, detectable by LIGO and by Advanced LIGO 
(O'Leary, O'Shaughnessy & Rasio 2007). Inspirals of a NS or of a stellar mass BH into an 
IMBH in GCs are detectable by Advanced LIGO, with rates up to tens per year (Brown 
et al. 2007; Mandel et al. 2008). Finally, ET might be able to observe GWs from IMBHs 
produced by seed BHs at high redshift (Sesana et al. 2009; Gair et al. 2009a, 2009b). 

In this Paper we extend the previous studies by Miller (2002) and by Mandel et al. (2008) 
to IMBHs hosted in YCs, making predictions for the detection rate of GWs from IMBHs in 
YCs. YCs are an ideal environment in which GWs from IMBHs can be studied, as YCs are a 
key place for the formation of IMBHs and because the density of YCs is expected to rapidly 
increase with redshift (Hopkins & Beacom 2006). For this reason, we pay particular attention 
to the integration with redshift of the detection rate. We consider, as representative of the 
science capability of first, second and third generation ground-based instruments Initial 
LIGO, Advanced LIGO and ET, respectively, and for space-borne instruments LISA. We 
note that space-based instruments with peak sensitivity in the frequency region ~ 0.1 Hz 
that fills the gap between LISA and ground based laser interferometers would be ideal for 
studying IMBHs, due to the emission frequency of systems in the mass range ~ 10 2 — 10 3 M Q . 
Such instruments, such as DECIGO (Kawamura et al. 2006), ALIA (Bender et al. 2005) 
and the Big-Bang-Observer (Phinney et al. 2003) are being investigated, but their timescale 
is currently very uncertain and for this reason we will not consider them in this Paper. 
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2. Results 

In this Section, we will focus on the properties of IMBHs in YCs and on their role as 
sources of GWs. In particular, we will consider two galaxies: the MW, for which we have 
the best observational data about YCs, and the Antennae, which host a large population of 
YCs. These two galaxies can also be considered as prototypical of, respectively, a bright disk 
galaxy with a typical star formation (SF) rate for the current cosmic epoch, and a merging 
system involving two massive disk galaxies undergoing enhanced SF, a configuration that 
is expected to be increasingly common at higher redshift. We will subsequently generalize 
our results and make predictions for the detection rate of GWs by LIGO, Advanced LIGO, 
LISA and ET. 



2.1. IMBH mergers in the MW 

First, we derive an estimate of the number of IMBHs which may exist in the YCs of 
the MW. As the YCs reside in the high-density regions of our Galaxy, their detection and 
observation are made difficult by Galactic absorption. Thus, only data for ~ 10 massive 
(> 1O 4 M ) Galactic YCs are available (Davies et al. 2007 and references therein). However, 
Gvaramadze, Gualandris & Portegies Zwart (2008) estimate that the total number of massive 
Galactic YCs is ~ 70 — 100, much higher than the number of the detected ones. The runaway 
growth of IMBHs is possible only in those YCs with half-mass relaxation time < 25 Myr 
(PZM02). Extrapolating from the properties of the YCs listed in table 1 of Portegies Zwart 
et al. (2002), we find that a fraction f tot ~ 0.5 — 1.0 of the Galactic massive YCs have 
th < 25 Myr and thus may host IMBHs. Therefore, ~ 30 — 100 IMBHs may exist at present 
in the YCs of the MW. In the following sections, we adopt the average value f to t = 0.75 as 
a fiducial value. 

IMBHs born via runaway collapse in the center of dense YCs likely form binaries with 
other stars or stellar mass BHs (see Colpi et al. 2003 and references therein). Such binaries 
harden due to three-body (or even four-body) interactions at a rate (Sigurdsson & Phinney 
1993; Colpi et al. 2003) 

o n -l in -6 -l ( m BH \ ( n c \ ( a \ ^Okms" 1 ^ 

u sb ~ 2 Ti G m BH n c a <r„ ~ 10 yr - 

d y V5OOM y V5xl0 5 pc- 3 y V0.4A.U./ \ (t c J 

(1) 

where G is the gravitational constant, is the mass of the IMBH, n c is the core density 
of the YC, a is the orbital separation of the binary and a c is the velocity dispersion in the 
core of the YC. In equation ([I]) we take as fiducial values for the MW n c = 5 x 10 5 stars pc~ 3 
and <7 C = 20 km s _1 , which are the average values for the YCs listed in table 1 of Portegies 
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Zwart et al. (2002). The masses of IMBHs formed by runaway collapse (mbh) likely range 
from 10 2 to 10 3 M Q (PZM02; Giirkan et al. 2006). Finally, a typical initial orbital separation 
of ~ 0.4 A.U. is assumed, according to Devecchi et al. (2007)1. 

Typically, a merger between the IMBH and its companion occurs after a few hundred 
three-body interactions (Miller 2002). Since it is hard to give a more accurate expression for 
the merger rate, in the following we will assume that the merger rate for a single IMBH binary 
is iVrg = 10~ 2 ^3b (Miller 2002), neglecting its possible dependence on various quantities, 
such as iiibh- Thus, the merger rate for a single IMBH binary in the Milky Way is v mTg ~ 
10" 8 (m BH /500 M ) (n /5 x 10 5 pc" 3 ) (a/0.4 A.U.) (20 km s -1 /^) yr" 1 . We expect a Galactic 
total merger rate ^ m rg,tot ~ 10~ 6 yr _1 , considering all the ~ 100 IMBHs hosted in the YCs 
of the MW. 



2.2. IMBH mergers in the Antennae 

The Antennae are a well-studied nearby galaxy pair which is undergoing merger. As 
the merger triggers the SF, the Antennae are richer in YCs than the MW. In fact, more 
than 1000 YC candidates have been observed in the Antennae (Mengel et al. 2005). The 
expected number of massive YCs in the Antennae can be derived with the same calculation 
used by Gvaramadze et al. (2008) for the MW. In particular, considering a SF rate (SFR) of 
7.1 M Q yr" 1 for the Antennae (Grimm, Gilfanov & Sunyaev 2003), assuming an upper limit 
of ~ 3 x 10 7 yr for the age of the YCs and using a power-law cluster initial mass function 
(MF) with slope 2 (Gvaramadze et al. 2008), we find that ~ 2100 massive YCs are expected 
to exist in the Antennae at present. Assuming that the half- mass relaxation time for the 
YCs in the Antennae is similar to that of the YCs in the M\vB, we expect that ~ 1000 — 2000 
IMBHs exist in the Antennae. Thus, the expected merger rate for IMBHs in the Antennae 



lr The fiducial value a = 0.4 A.U., adopted in our calculations, is the typical orbital separation of an 
IMBH— BH (or IMBH— NS) binary which is already hard (i.e. whose binding energy is larger than the 
average kinetic energy of a star in the cluster, see e.g. Heggie 1975), but which is still in the regime where 
the hardening due to gravitational encounters is more important than that due to gravitational wave emission. 
For more details about the derivation of a, see Table 1, equations 1 — 3 and the Appendix A of Devecchi et al. 
(2007). We note that a more realistic approach should consider a distribution of values for a, depending on 
tobhj on the mass of the companion, on n c , on a c , and on the history of previous gravitational encounters. 
Here, we simply adopt the peak of the distribution of a for mBH ~ 300 M Q , as derived by Devecchi et al. 
(2007), in order to provide an approximate estimate of v^. 

2 This assumption is reasonable, as various properties of the YCs in the Antennae (e.g. velocity dispersion, 
total cluster mass, etc.) are similar to the ones of the YCs in the MW (Mengel et al. 2002; de Grijs et al. 
2005). 
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is z/ mrgitot « 10 yr , i.e. a factor of > 10 higher than for the MW. 



2.3. GWs from IMBHs in the local Universe 

From the calculations reported in the previous Sections, we can predict a detection 
rate of GWs emitted by merging IMBHs in the local Universe, expected for LISA, LIGO, 
Advanced LIGO and ET. In particular, we consider all the mergers involving IMBH— BH and 
IMBH— NS binaries on the past light cone, within a sphere of fixed radius, derived from the 
detection range of a given instrument, through the expression of the maximum observable 
redshift z max , see the description below. Under the assumption that all the IMBHs have the 
same mass, such detection rate can be expressed as (Miller 2002): 



3 Zmax(m B H,mco) / Z \ 2 

\0 / 

where c is the light speed, H (= 72 km/s Mpc -1 , Spergel et al. 2007) the Hubble constant, 
/tot ~ 0.75 is the fraction of massive YCs hosting an IMBH (see Section 2.1), v mTg (mBii) is 
the merger rate per IMBH (see Section 2.1), E(z) = [(1 + z) 3 n M + fi A ] 1/2 (where fi M = 0.27 
and f2yv = 0.73, Spergel et al. 2007). The factor ^ = (1 + z)~ l accounts for the difference 
between the time in the rest frame of the source (t c ) and of the observer (t ). 

In equation (j2J), nyc{z , rriBn) is the comoving number density of YCs which are suffi- 
ciently massive to host an IMBH. Assuming that a constant fraction of stars (/sfc) forms 
in clusters, nyc (^> ^bh) can be approximately estimated as: 



n Y c{z,m Bn ) = 4.5 x 10~ 2 Mpc 3 



P*0) 



1.5 x lO^Moyr^Mpc" 
-i 



( t max \ / /surv \ / (m Y c) \ f /SFC \ ( ZyC^BhA , , 

X \W 9 yr) V 10 " 2 / V260M.J V 0.8 J \ 10~ 4 )' [ ) 

where t max is the maximum lifetime of a YC: we adopt t max = 10 9 yr, as we account for the 
fact that the IMBH— BH and IMBH— NS mergers occur ~ 2 x 10 8 yr after the formation of 
the central binary. We must also include a correction f SUTV which represents the fraction of 
YCs which avoid disruption and survive up to t max = 10 9 yr. / surv is quite uncertain; we 
adopt a conservative value / surv = 10 -2 from the literature (Lada & Lada 2003; Fall et al. 
2005). For a discussion about the possible mechanisms of cluster infant mortality, see, e.g., 
Gieles (2010) and references therein. 
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p # (z) = 10^ lo §io i s the comoving density of SFR (where a = -1.82 and (3 = 3.28 

for z < 1.04, and a = -0.724 and /3 = -0.26 for 1.04 < z < 4.48, Hopkins & Beacom 2006). 
These values of p*(z) are based on the most updated availlble datg 

/sfc is the fraction of 

SF which occurs in clusters (/sfc — 0.7 — 0.9, Lada & Lada 2003). We assume that /sfc 
is constant with redshift, although this cannot be proved on the basis of currently available 
data. 

(myc) ~ 260 M Q is the average mass of a YC, and has been derived assuming that the 
masses of the YCs range from m Y c,min ~ 20 M to m Y c,max ~ 1O 7 M and are distributed 
according to a MF 4^ oc m~ 2 (Lada & Lada 2003). Finally, /yc^bh) is the fraction of YCs 
which are massive enough to host IMBHs: 



/yc depends on jtibh, as the simulations show that the IMBH mass scales with the cluster 
mass (ttt-bh ~ 10 _3 mYc ; PZM02). In the following calculations, we assume that myc, mm = 
20 M , that m Y c,max = 1O 7 M and that ^ oc m~ 2 (Lada & Lada 2003). Under these 
assumptions, /yc ( m BH) g° es from ~ 2 x 10~ 4 (for an IMBH mass m B H = 100 M ) to 
~ 1.8 x 10~ 5 (for an IMBH mass m BH = 1000 M ). 

In equation ([2]), -z max (^BH ; m co) is the maximum redshift at which an event can be 
detected with a sky-location and orientation averaged signal-to-noise ratio (SNR) > 10 by 
a single interferometer. In observations with a network of instruments, the signal-to-noise 
ratio scales as the square root of the number of instruments, and in this respect the results 
presented here should be considered as conservative. The maximum redshift depends on the 
mass of the IMBH meH and on the mass of the companion which merges with the IMBH 
(m co ), as well as on the sensitivity of the instrument. The expressions adopted to calculate 
2 max( m BHi m co) (for LIGO, Advanced LIGO, LISA and ET) are summarized in Appendix A. 
As a model for the gravitational waveform, we consider the analytical phenomenological 
inspiral-merger-ringdown waveform approximant for non- spinning BHs in circular orbits 
derived by Ajith et at. (2008a). This waveform model together with the Effective- One- 
Body- Numerical- Relativity (EOBNR) waveform family (Buonanno et al 2007) provides a 
prescription for the signal from the whole coalescence of binary systems. These waveform 



3 We note that the assumption that the YC density is proportional to the SFR is quite inconsistent with 
the fact that we assume a lifetime t max = 10 9 yr for the YCs, likely overestimating high-redshift clusters 
(z ~ 4) by < 30 per cent. On the other hand, it is quite hard to refine this model (given our poor knowledge 
of high-redshift YCs). Furthermore, values of the redshift higher than 1 are important only for ET (see 
FigfT]). Finally, we adopted a conservative value of / sm v (= 10~ 2 ), in order to avoid overestimating nyc- 
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approximants have both been calibrated on full numerical relativity simulations for binaries 
for approximately equal mass systems, and yield signal-to-noise ratios that are consistent 
within w 25%. However, in this paper we will apply the waveforms on a much larger mass- 
ratio regime that extends up to ~ 10~ 2 . For this mass-ratio we still do not have reliable 
inspiral- merger-ring-down waveforms, and the phenomenological and EOBNR approximants 
yield significantly different signal-to- noise ratios, that differ by a factor « 2—10 depending 
on the mass ratio, with the phenomenological waveforms producing the larger signal-to-noise 
ratio, see e.g. Fig. 2 of Gair et al (2009b). As the detection rate scales as the cube of the 
maximum distance at which a source can be detected, the amplitude differences lead to rate 
uncertainties of a factor 10 or larger. Furthermore if the BHs are (rapidly) spinning, the 
signal-to- noise ratio at which a source is observed can be significantly (by a factor of ~ 2, 
see e.g. Ajith et al. 2009) affected. However, the lack of astrophysical predictions for the 
likely spins of IMBH binaries and the lack of full coalescence waveform approximants for 
generic spin magnitudes and geometries (see however Ajith et al. 2009; Pan et al. 2009) 
prevent us from considering this possibly important physical effect. The resulting values of 
2 m ax(^BH, ^co) are shown in the upper panels of Fig. [TJ In this Fig. we consider mergers of 
the IMBHs with two different compact objects, i.e. NSs (left-hand panel) and stellar mass 
BHs (right-hand panel)@. It is evident that there is a large difference between ET and the 
other interferometers. In fact, while for ET Zmax^BH, ^co) is always higher than 0.5 and 1 in 
the case of NSs and stellar mass BHs, respectively, for Advanced LIGO Zmax^BH, ^co) < 0.2 
and < 0.5 in the case of NSs and stellar mass BHs, respectively. For the other interferometers 
z mffl(( m BH, m co) is even smaller. 

The detection rates R derived from equation ([2]) for events detectable with Advanced 
LIGO, LISA and ET are shown in the bottom panels of Fig. [IJ In particular, the detection 
rates R for mergers involving IMBH— NS and IMBH— BH binaries are shown in the left-hand 
panel and in the right-hand panel, respectively. The different lines show the values of R for 
different detectors, as a function of the IMBH mass, ranging from 10 2 to 10 3 M . The results 
obtained for LIGO are not shown in Fig. [TJ as they are orders of magnitude lower than those 



In our calculations we assume a mass of m co = 1.4 and m co = 10 M for the NSs and the stellar mass 
BHs, respectively. Furthermore, we assume that all the IMBHs formed in YCs merge with a compact object. 
This assumption is justified by the strong mass segregation occurring during the runaway collapse (PZM02). 
The strong mass segregation and the concentration of the most massive stellar objects at the center of the 
cluster are confirmed by the observed mass distribution in the core of some Galactic globular clusters (where 
the mass distribution within the core is 4^ oc m a , with a = 3 — 8, see e.g. Prior et al. 1986 and Monkman 
et al. 2006). These evidences, combined with the fact that hard binaries (i.e. binaries with binding energy 
larger than the average kinetic energy of a cluster star) tend to exchange, during three-body encounters, with 
the most massive possible companion (see e.g. Sigurdsson & Phinncy 1995), suggest that most of IMBHs 
formed in YCs merge with compact objects. 
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obtained for the other interferometers. In particular, the rate of IMBH— NS (IMBH— BH) 
mergers detectable with LIGO is lower than 5 x 10 -5 (1 x 10~ 3 ) per year, even in the most 
favorable case (corresponding to m^n = 100 M Q ). 

Fig. [1] shows that ET (dotted line, red on the web) is expected to detect a large number 
of events per year: R > 10 yr _1 and R > 60 yr _1 in the case of IMBH— NS and IMBH— BH 
mergers, respectively. Thus, ET will be essential to assess the existence of IMBHs in the 
Local Universe, and to distinguish among different theories of IMBH formation. Advanced 
LIGO (solid black line) is expected to detect a few mergers involving IMBH— BH binaries 
per year and less than one merger involving IMBH— NS binaries per year. Advanced LIGO 
is particularly suitable for observations of low-mass IMBHs (< 300 M ) merging with stellar 
mass BHs. In the case of m B H = 100 M Q , Advanced LIGO is expected to detect up to ~ 5 
IMBH— BH mergers per year. 

On the other hand, LISA (dashed line, green on the web) is more efficient in detect- 
ing massive (> 500 M Q ) IMBHs merging with stellar mass BHs, due to its sensitivity to 
frequencies lower than 1 Hz. However, LISA is not particularly suited for detecting GWs 
from IMBHs, as less than one event per year is expected to be detected, even in the most 
optimistic case (i.e. for ttibh = 1O 3 M ). In the case of IMBH— NS mergers, LISA is very 
unlikely to observe any such system during the mission lifetime (nominally set to 5 years), 
even for ttibh = 10 3 M Q ; the detection rate is in fact < 10 _4 yr _1 . 

In equation §2§ we assume that all the IMBHs have the same mass. However, this 
is an over-simplification, as IMBH masses will be distributed according to a certain MF. 
Accounting for such MF, equation ([2]) becomes: 

m 2 

R = 4 Ti ( — ) ftnt I / - dm BH I 




in-2 

(wbh) 

»'l 

_dN_ 

im B ] 

IMBH mass, respectively. 



dm B H 

mi 

ZmaxfmBE, m C o) 




n Y c(z,m BH ) , 
(l + z)E(z) Z 



dN A (rt 

■ dm B H, (5) 



dm BH 



where is the MF of the IMBHs, whereas rrii an d nri2 are the minimum and the maximum 



The MF of IMBHs is unknown. However, as we already mentioned for the estimate of 
/yc(^bh), PZM02 show that there is a correlation between the mass of the IMBH and the 
mass of the host YC (myc), and in particular that the mass of the IMBH is wbh ~ 10 _3 mYc- 
Since the YCs in the MW have a MF 4^- oc m~ 2 (Lada & Lada 2003), we can assume 
that the IMBHs have the same MF. Adopting this MF and assuming mi = 100 M and 
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m 2 = 1000 M Q , we find (m BH ) = 256 M . 

Integrating equation (jSJ) for these values of the MF and of mi and m 2 , we obtain the 
detection rates R shown in Table 1. In Fig. [JJwe also plot the results of equation (j5J), using 
a filled square (red on the web), a filled black circle and a filled triangle (green on the web) 
for ET, Advanced LIGO and LISA, respectively. These points are drawn, for simplicity, in 
correspondence of ttt-bh = 256 M , but they have been obtained integrating equation over 
a Lada & Lada (2003) MF, with an average IMBH mass (771bh) = 256 M . We note that 
the results obtained from equation (|S]), under these assumptions, are very similar to those 
obtained from equation (J2J), assuming that all the IMBHs have mass ttibh = 256 M Q . 

Fig. [1] and Table 1 show that, if the IMBHs are distributed according to a Lada & 
Lada (2003) MF, ET is expected to detect a large number of events per year, involving both 
IMBH-NS (R ~ 200 jt' 1 ) and IMBH-BH (R ~ 400 yr" 1 ) binaries. Advanced LIGO is 
expected to detect ~ 4 mergers involving IMBH— BH binaries and ~ 0.2 mergers involving 
IMBH— NS binaries per year. LISA is not suited for detecting IMBHs distributed according 
to the Lada & Lada (2003) MF, as such MF favors relatively 'light' IMBHs. 

In Appendix B we derive an alternative calculation of R, based on equation (14) of Gair 
et al. (2009b). In Fig. [TJ the results of the alternative calculation reported in Appendix B are 
shown as open points: an open square in the case of ET, an open circle for Advanced LIGO, 
and an open triangle for LISA. The results of the alternative calculation are similar (within 
a factor of 5) to those of equation (J3J. Given the large uncertainties in both calculations, 
the results are consistent between them. Furthermore, a more simplified, order-of-magnitude 
derivation of the detection rate R is provided in Appendix C. We stress that the calculation 
in Appendix C is much more approximated than the one presented in this Section and the one 
reported in Appendix B, but it is a useful tool, in order to estimate the order of magnitude 
of R with 'back of the envelope' calculations. 

Our results are consistent with those indicated by previous studies. In particular, in the 
case of LISA, our estimates of the detection rate R are approximately one order of magnitude 
higher than those obtained in Will (2004), which adopts the noise curves by Larson et al. 
(2000), and a factor of < 5 lower than those derived by Miller (2002), based on the noise 
curves in Flanagan & Hughes (1998a, 1998b). The results obtained for Advanced LIGO are 
consistent with the recent study by Mandel et al. (2008). In the case of ET, as we discuss 
in Appendix C, the results of our approximate calculation in equation (ICip are very similar 
(within a factor of two) to those reported in Table 1 of Gair et al. (2009b). The results 
derived from the more refined equation (jSJ) and from the alternative model in equation (IB4|) 
are a factor of 4 — 10 larger than those reported in Gair et al. (2009b), mainly due to the 
fact that we integrate over the IMBH mass function and that we account for the dependence 
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of iiyc on the redshift. 

In conclusion, GWs from IMBHs are hardly detectable with LIGO and LISA. Instead, 
Advanced LIGO has chances of detecting GWs from IMBHs in clusters. ET is expected to 
observe hundreds of IMBHs formed via the runaway collapse mechanism. 

3. Summary 

In this Paper we study the occurrence of mergers between IMBHs and compact objects 
(NSs and stellar mass BHs) in YCs. These are found to be important sources of GWs. 
Our study shows that GWs from IMBHs are unlikely to be detected with first generation 
instruments, such as Initial LIGO, and even with LISA. Advanced LIGO offers instead the 
ability of observing these sources: a few merger events are expected to be detected by 
Advanced LIGO in 1-year integration. ET is far the best interferometer projected up to now 
to study GWs from IMBHs: ~ 10 to ~ 300 (~ 60 to ~ 600) merger events of IMBH-NS 
(IMBH— BH) binaries might be detected in 1-year integration, according to the scenario of 
runaway collapse for IMBH formation. Thus, ET will be a powerful tool to check the runaway 
collapse and the other scenarios for IMBH formation. However, our estimates are affected 
by large uncertainties (more than a factor of 10), because of our insufficient knowledge of 
masses and density of IMBHs in the local Universe. In particular, the fraction of IMBHs per 
cluster / to t is highly uncertain, as it is based only on the properties of the ~ 5 Galactic YCs 
for which enough data are available. Furthermore, / to t = 0.75 adopted in this Paper must 
be considered as an upper limit, as we are assuming that all the YCs which can undergo 
runaway collapse host an IMBH, and that all the mass involved in the collapse ends up 
into the IMBH. There are no evidences that the runaway collapse necessarily leads to the 
formation of an IMBH (PZM02) and there are hints that a large fraction of the mass is lost 
due to winds and/or to recoil (Gaburov, Lombardi, Portegies Zwart 2010). Thus, there is no 
evidence that the MF of IMBHs is the same as that of the host YCs. Furthermore, even the 
density (wyc(z, tibh)) and the MF of YCs are uncertain: in order to derive such quantities, 
we extrapolated to massive (> 1O 4 M ) YCs various relationships (e.g. the slope of the MF) 
derived by Lada & Lada (2003) for a sample of smaller clusters (20 — 1000 M ). 

Finally, the models adopted to derive the instrumental range, and the corresponding 
z max , of the considered interferometers assume that the binary has zero eccentricity, the 
BHs are non-spinning and consider waveforms in a mass-ratio regime in which we still do 
not have reliable signal approximants. These are all factors that may significantly alter the 
results presented in this paper, in particular should YCs produce rapidly spinning and/or 
highly eccentric binaries. Full inspiral-merger-ringdown gravitational-waveform families for 
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(generic) binary systems with a mass ratio of the order of 1— to— 100 are not available at 
present and we will be able to address this problem more rigorously only in the future. 

We also note that, although we have consistently taken into account the non-linear 
dynamics of the binary during the final merger by adopting full inspiral-merger-ringdown 
waveforms (Ajith et al 2008a), calibrated on actual numerical relativity simulations of the 
coalescence of BH binaries, we have however neglected in the merger rate estimates the 
possibility that the GW recoil (see, e.g. Baker et al. 2008; Lousto et al. 2010) may impart 
to the binary a sufficient recoil velocity to eject it from the cluster. This would prevent 
future mergers, and in this respect, the results presented in this paper should be regarded as 
upper-limits, as GW recoil can only reduce the merger rate. On the other hand, assuming 
a typical escape velocity of 20 km s _1 from the core of the host cluster (see e.g. Colpi et 
al. 2003), the GW recoil (calculated according to Baker et al. 2008) should not be able to 
eject IMBHs more massive than ~ 250 M (~ 35 M ) merging with a 10 M Q stellar mass 
BH (with a 1.4 M NS). We notice, however, that the available expressions of the recoil 
velocity were obtained from simulations with mass ratio m co /mBH = 1 — 1/3 (e.g Baker et 
al. 2008), whereas we consider systems for which m co /m B H = 0.1 — 0.0014. No simulations 
have been carried out for such extreme mass ratio. Furthermore, even in a pessimistic 
scenario, accounting for GW recoils should reduce the detection rate R by at most a factor 
of a few, as in our model most of IMBHs undergo < 5 merger events within the lifetime of 
the cluster. 

All these caveats must be taken into account, when comparing the rates R derived in 
this Paper with the forthcoming observational data. Yet our results will provide the basis for 
future studies exploring these issues further with improved data on the nature and frequency 
of YCs. 

We thank the referee for his/her critical reading of the manuscript. MM acknowledges 
support from the Swiss National Science Foundation, project number 200020-109581/1. 
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A. Method to estimate z max 

We define z ma x(^BH; m co) as the maximum redshift at which an event can be detected 
with a sky location and orientation averaged signal-to-noise ratio (SNR) > 10 by a single 
interferometer. The coherent SNR at which a gravitational wave h(t) can be detected by an 
instrument characterized by a (zero-mean) Gaussian, stationary noise with a one-sided noise 
power spectral density S n (f), is given by 

SNR 2 = (h\h). (Al) 

Here, (.|.) is the noise- weighted inner product, defined as (Cutler and Flanagan 1994) 

{a]b) ^f mm+mhfi df , (A2) 

where a(f) and b(f) are two generic functions defined in the Fourier domain. 

For the waveform h(t) we use the phenomenological inspiral-merger-ringdown waveform 
model by Ajith et al. (2008a) for non-spinning BHs in circular orbit. For an optimally 
orientated binary the waveform emitted by during the whole coalescence is described by: 

u(f) = A eS (f) e™^, (A3) 



This preprint was prepared with the AAS IATgX macros v5.2. 
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where 

{(// /merg) if / ^ /merg 

(///merg)^ 2/3 if /merg < / < /ring (A4) 
W£(/, /ring, 0-) if /ring < / < /cut 

Expressions for / mcrg , / ring , / cut , (/), C, 10 and £(/, / ring , a) are given in equations (4.14- 
4.19) and Tables I-II of Ajith et al. (2008a). The fitting coefficients reported in Ajith et al. 
(2008a) have been revised in Ajith (2008b), but have no effect on the actual signal-to-noise 
ratio, and we have therefore adopted the original values. The angle-averaged signal-to-noise 
ratio is obtained by dividing the SNR, equation OAip from an optimally oriented source, 
described by equation (IA4[) . by a factor 2.26, which accounts for the varying response of GW 
instruments to sources in different locations of the sky and with different orientations of the 
orbital angular momentum: 

(SNR) = M {M) 

In observations carried out with a network the total coherent network SNR is (SNR) 2 = 
^fc(SNR) 2 ,, where (SNR) 2 is the signal-to-noise ratio at each instrument. We caution the 
reader that if sources are not uniformly distributed in distance, then considering angle- 
averaged signal-to-noise ratios (equation (1A5|) ) introduces some errors in the estimate of the 
detection rates; however, considering the other (large) uncertainties (from the astrophysics 
and waveform modelling) that enter the computation of this quantity, the simplification 
introduced by equation flA5j) has a negligible effect on the final results. 

Using these equations, we can derive the luminosity distance -Di(2 ma x("iBH, ^co)) at 
which an event can be detected with a SNR > 10 by a certain interferometer. We can, thus, 
derive z m ax( m BH, m co) by inverting the expression of the luminosity distance in the A Cold 
Dark Matter (CDM) model: 

Ac^max^BH, ma,)) = tH 1 + 2 max(™BH, ^co)) / -=^r (A6) 

Ho Jo E(z) 



The range of an instrument, equation fl A6 j) and therefore the IMBH binary detection 
rates are entirely determined by the noise spectral density S n (f), see equation (I AID . For 
ground-based interferometers we consider representative sensitivity curves for the three gen- 
erations of instruments. For instruments now in operation (first generation) we adopt the 
initial LIGO design sensitivity curve, that well approximates the sensitivity achieved during 
the last science run (Abbott et al. 2009). Consistently we set the low frequency cut-off, the 
minimum of integration in equation flA2j) . to f mm = 40 Hz; the upper frequency cut-off is 
irrelevant, as the ring-down signal for IMBH mass-range of interest is de facto zero in the 
high-frequency region of the instrumental sensitivity window / ~ 1 kHz. 



- 18 - 



For second-generation (or advanced) interferometers we use the broad band target design 
sensitivity curve of Advanced LIGO with a low frequency cut-off / m ; n = 10 Hz (Adhikari et 
al, 2006); Advanced Virgo is expected to operate on the same timescale and has similar noise 
performance. 

For third generation interferometers, that are currently undergoing conceptual design 
studies, we adopt the noise curve of a single right-angle ET instrument and a low frequency 
cut-off / min = 1 Hz (Hild et al. 2008). We note that other configurations, such as the 
Xylophone (Hild et al. 2010), have been proposed and are currently under study. 

For LISA we adopt the current best estimate of the instrumental noise spectral den- 
sity, see Barack and Cutler (2004). As we show below, the frequency range relevant for 
observations of IMBH binaries is always above several mHz; as a consequence, the confusion 
noise generated by stochastic foregrounds of close white-dwarf binaries (and possibly extreme 
mass-ratio inspirals) does not contribute to the total noise budget and we therefore ignore 
it in the calculation. In the case of low-frequency observations, the binary system lifetime 
may be longer than the mission lifetime. In fact, the time to coalescence (at the leading 
quadrupole Newtonian order) for a binary radiating at frequency / is 

r(/) = 5c 5 G- 5 / 3 (8 7r/)- 8 / 3 /i- 1 m-f , (A7) 

where G is the gravitational constant, c is the speed of light, and m tot = m B H + m co and 
\i = mBH m co/ m tot are the total and reduced mass, respectively, of the coalescing binary. 
Moreover, LISA will not be able to observe the final merger-ringdown phase, as it takes 
place at frequency higher than the observable window. Here we assume that the highest 
frequency that LISA can observe is 1 Hz and the last stable orbit around a Schwarzschild 
BH corresponds to the GW frequency ~ 4 (10 3 M Q /m tot ) Hz. The integration limits in 
equation flAip are obtained with the following procedure. The time of observability T obs for 
a merger event is given by (Peters 1964) 

Tobs = r(f min ) - r(/ max ) . (A8) 

From equation (IA8I) . we can derive / min as a function of / max , assuming T obs = 5 yr. Fig. [2] 
shows / m i n as a function of /max- We note that / m in is almost constant for / max >3x 10~ 2 
Hz. We then calculate (SNR) from equation (1A5I) for different values of the couple /min, 
/max- Fig. [3] shows (SNR) as a function of the couple / m in, /max- We find that (SNR) is 
maximum when we take / max = 1 Hz and we consistently derive / min from equation (1A8|) . 
The results presented in Table 1 and in Fig. [TJ are obtained with such choice of f max = 1 Hz 
and / m i n derived from equation flA8j) . As a consequence the results presented for LISA in 
this Paper should be considered as upper limits, assuming that the binary is in the optimal 
stage of the merger for detection. However, from Fig. [3] we also note that (SNR) is almost 
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constant for / max > 10~ 2 Hz. Thus, we expect that the upper limit for LISA is not far from 
the average value. The integration frequency range also justifies neglecting the confusion 
noise from unresolved stochastic foregrounds. 



B. Alternative derivation of R 

Gair et al. (2009b) report in their equation (14) a derivation of the detectable merger 
rate for IMBH— IMBH binaries. In the current paper, we do not consider the mergers of 
IMBH— IMBH binaries, because the possibility of forming more than one IMBH in the same 
cluster is still debated (see e.g. Giirkan, Fregeau & Rasio 2006). However, it is possible to 
perform a calculation similar to equation (14) of Gair et al. (2009b) for IMBH— BH and 
IMBH— NS binaries, as follows. 

*=/ *"B„/ d ^. r 8 ^ ^^. (BH 



mi J 



dm B H dt c dV c dt Q dz 



where mi and ni2 are the minimum and the maximum IMBH mass, respectively; ^ = 
(1 + z)^ 1 (since t e and t Q are the time in the rest frame of the source and of the observer, 
respectively); V c is the comoving volume; A mcrg is the number of mergers. Assuming that 
m BH = 10- 3 m YC (PZM02), 



d 3 A merg _ 3 d 3 A YC 

— — ZVrgl^BHj Jtot £max Jsurv -W , [dZ) 

dmBH dt c d V c dm yc dr e d V c 



where N YC is the number of YCs. Following equation (12) of Gair et al. (2009b) 

P*(z)my 2 c , (B3) 



d 3 A YC fi'Yc(^Yc) 



dm YC dt e dV c In (m YC ,max/ m YC ,, 

where p*(z) is the comoving density of SFR (see Section 2.3) and 5'yc( ?ti yc) is the fraction 
of the total stellar mass which is formed in YCs of mass m Y c- 

Substituting equation ( 1B3|) into equation f ]Blj) . and assuming that ^mrg(^BH) = 10~ 2 ^(mBH") 
(see Section 2.1) and mBH = 10 _3 m Y c, we finally obtain: 

2 7T G /tot imax /surv , n _5 -1 /" mYC ' maxBH dm YC 

R = r~< / v 10 Uc a a c 9yc / 

m (m YC .max /m YC m YC 

Jo (1 + z) dz' 

where m Y c, m inBH and m Y c, m axBH are the minimum and the maximum YC mass in order to 
form an IMBH in the considered mass range, respectively (we adopt m Y c, m inBH = 10 5 M and 
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TtYC,maxBH = 10 6 M , corresponding to meH between 10 2 and 10 3 M Q ). As in Section 2.3, 
m YC,mm = 20 M Q and 77iYC,max = 10 7 M are the minimum and the maximum YC mass. 
In equation (IB4I) we assume that n c , a, a c and gyc do not depend on the cluster mass. 
Adopting n c = 5 x 10 5 pc~ 3 , a = 0.4 AU, a c = 20 km s _1 and gyc = 0.8 (the same values as 
in the alternative calculation reported in the main text), we derive the values of R shown in 
Fig. [1] as open points (open triangles for LISA, open circles for Advanced LIGO and open 
squares for ET). The difference between open and filled points (derived from equation IB4I 
and from equation respectively) is within the large uncertainties of the two calculations. 
In Table 2, we report the values of R derived from equation (IB4j) . for a comparison with 
those listed in Table 1: differences are less than a factor of 7 in the case of LISA and less 
than a factor of 2 for the other instruments. 



C. Approximate derivation of R 

An order-of-magnitude estimate of R may be obtained with a much simpler calculation 
than those reported in Section 2.3 and in Appendix B. In particular, assuming that all the 
IMBHs have the same mass wibh, we can write: 

fi^/, o , " Y ° (2 l 'T ) ^). < ci > 

i mrg ^m B Hj 

where nyc{z = 0,?tibh) is the comoving density of YCs which are sufficiently massive to 
host an IMBH, defined by the equation Q, calculated at z = 0. We note that we are 
considering a lower limit of the YC density, as the value of nyc(z, mbh) is minimum for 
z = 0. T mrg = z/"J g = 2 x 10 8 yr (256M /m BH ) (5 x 10 5 pc" 3 /n c ) (<r c /20 kms" 1 ) is the 
timescale for a merger. Finally, V c (z max ) is the comoving volume up to redshift z ma x(^BH) 
(with z max (mBn) derived as described in Appendix A). The results of equation fIClj) as a 
function of the fiducial IMBH mass are shown in Fig. HI Table 3 reports the value of R 
for a fixed IMBH mass mbh = 256 M Q . We note that the detection rates predicted by 
equation ( IClj) are a factor of ~ 2 — 10 lower than those derived with the other two, more 
refined, methods (see Fig. [Q. This is likely due to the fact that in equation fICip n YC does 
not increase with redshift. For this reason, equation ( 1C1|) is a sort of lower limit for our 
estimates. 

Finally, Table 4 is focused on the case of ET, and shows both the quantities adopted 
in equation (ICip and the derived detection rates for some choices of m B H an d m co . Table 4 
was introduced for a comparison between our procedure and the similar one reported in 
Section 3.1 and in Table 1 of Gair et al. (2009b). We note that the values of the comoving 
density of young clusters which are sufficiently massive to host an IMBH (nyc), adopted in 
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our calculations, are a factor of 3 — 40 lower than those adopted by Gair et al. (2009b), who 
consider the comoving space density of globular clusters and assume that ~ 10 per cent of 
them host an IMBH. On the other hand, the usage of Ajith et at. (2008a) model for the 
gravitational waveform leads to an overestimate of z max by a factor of 2 — 3 with respect to 
the EOBNR models adopted by Gair et al. (2009b, see the discussion in the main text for 
details). For these reasons, our estimates of R are very similar (within a factor of two) to 
those reported in Table 1 of Gair et al. (2009b). 
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Fig. 1. — Top panels: z max (mBK, wico) a s a function of the IMBH mass ttt,bh for IMBHs 
merging with NSs (left-hand panel) and with stellar mass BHs (right-hand panels). Bottom 
panels: detection rate R as a function of the IMBH mass mbh for IMBHs merging with 
NSs (left-hand panel) and with stellar mass BHs (right-hand panels). For all the panels, 
dotted line (red on the web): events detectable by ET; solid black line: events detectable 
by Advanced LIGO; dashed line (green on the web): events detectable by LISA. The events 
detectable by LIGO are not shown in this Fig., because they are more than one order of 
magnitude below the limits of the y— axes. For the same reason, Zmax^BH, ^co) and R for 
IMBH— NS mergers are not shown in the case of LISA. In the bottom panels, the filled square 
(red on the web), the black filled circle and the filled triangle (green on the web) are the 
detection rate R in the case of ET, Advanced LIGO and LISA, respectively, obtained from 
equation (jSJ), assuming a Lada & Lada (2003) MF for the IMBHs. These points are drawn, 
for simplicity, in correspondence of tbbh = 256 M , but they have been obtained integrating 
equation (jSJ) over a Lada & Lada (2003) MF, with an average IMBH mass (j71bh) — 256 M 
(see Section 2.3 for details). In the bottom panels, the open square (red on the web), the 
black open circle and the open triangle (green on the web) are the detection rate R in the 
case of ET, Advanced LIGO and LISA, respectively, obtained from equation (1B4|) . assuming 
a Lada & Lada (2003) MF for the IMBHs. 
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Fig. 2. — / min as a function of / maX ; i n the case of LISA, for a binary at z — 0.1 with 
m co = 10 M Q . Solid black line: wbh = 100 M ; dotted line (red on the web): ttibh = 
200 M Q ; dashed line (green on the web): mbh = 500 M ; dot-dashed line (blue on the web): 
m BH = 1000 M . 
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Fig. 3. — (SNR) as a function of f mSuX , in the case of LISA, for a binary at z — 0.1 with 
m co = 10 M Q . Solid black line: wbh = 100 M ; dotted line (red on the web): mbh = 
200 M ; dashed line (green on the web): mbh = 500 M ; dot-dashed line (blue on the web): 
m BH = 1000 M . 
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Fig. 4.- Left-hand (Right-hand) panel: detection rate R for IMBH-NS (IMBH-BH) 
mergers derived with the order-of-magnitude calculation reported in equation ( IClj) . Dotted 
line (red on the web): events detectable by ET; solid black line: events detectable by Ad- 
vanced LIGO; dashed line (green on the web): events detectable by LISA. For comparison, 
we include also the results of equation (J5]). In particular, the filled square (red on the web), 
the black filled circle and the filled triangle (green on the web) are the same as in Fig. [H i.e. 
the detection rate R in the case of ET, Advanced LIGO and LISA, respectively, obtained 
from equation (jSJ), assuming a Lada & Lada (2003) MF for the IMBHs. 
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Table 1: Detection rates from equation (jSJ) assuming a Lada & Lada (2003) MF and /tot = 
0.75 for the IMBHs (see Section 2.1 and 2.3 for details). 





LIGO 


Advanced LIGO 


LISA 


ET 




R [yr- 1 ] 


R [yr- 1 ] 


R [yr- 1 ] 


R [yr- 1 ] 


NS 


2 x 10" 5 


0.2 


2 x 10" 5 


200 


stellar mass BHs 


3 x 10" 4 


4 


0.01 


400 
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Table 2: Detection rates from equation (IB4I) assuming a Lada & Lada (2003) MF and f to t = 
0.75 for the IMBHs (see Appendix B for details). 





LIGO 


Advanced LIGO 


LISA 


ET 




R [yr- 1 ] 


R [yr- 1 ] 


R [yr- 1 ] 


R [yr- 1 ] 


NS 


2 x 10" 5 


0.3 


1 x 10" 4 


300 


stellar mass BHs 


3 x 10" 4 


6 


0.07 


750 
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Table 3: Detection rates from equation (IC1I) assuming wibh = 256 M Q and f to t = 0.75 for 
the IMBHs (see Appendix C for details). 





LIGO 


Advanced LIGO 


LISA 


ET 




R [yi- 1 ] 


R [yr- 1 ] 


R [yr- 1 ] 


R [yr- 1 ] 


NS 


2 x icr 6 


0.2 


1 x 10" 5 


30 


stellar mass BHs 


3 x 10" 5 


2 


7 x 10" 3 


100 
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Table 4: Quantities used in equation (ICip and detection rates R, derived from the same 

equation, in the case of ET (see Appendix C for details). 

m BH [Mp] w co [M Q ] z max n YC [Mpc~ 3 ] T mrg [yr] V c [Mpc 3 ] g [yr' 1 ] 

100 1.4 2.1 0.09 5.2 x 10 s 6.1 x 10 11 80 

100 10 4.4 0.09 5.2 x 10 8 1.7 x 10 12 200 

300 1.4 1.2 0.03 1.7 x 10 8 2.1 x 10 11 30 

300 10 2.3 0.03 1.7 x 10 8 6.9 x 10 11 90 

1000 1.4 0.5 0.008 5.2 x 10 7 2.9 x 10 10 3 

1000 10 1.0 0.008 5.2 x 10 7 1.5 x 10 11 20 



